NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 

No. 1818 


EFFECT OF AUTOMATIC STABILIZATION ON THE LATERAL 

OSCILLATORY STABILITY OF A HYPOTHETICAL 

AIRPLANE AT SUPERSONIC SPEEDS 

By Leonard Sternfield 

Langley Aeronautical Laboratory 
Langley Air Force Base, Va. 



EEBhhTD 



TECH LIBRARY KAFB, NM 


91 

991! 

911 

9EU 

111 



□: 

LMMI 

33 

3 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE NO. l8l8 

EFFECT OF AUTOMATIC STABILIZATION ON THE LATERAL 
OSCIILATORY STABILITY OF A HYPOTHETICAL 
AIRPLANE AT SUPERSONIC SPEEDS 
By Leonard Sternfield 

SUMMARY 


A theoretical investigation has "been made to determine the effect 
of automatic stabilization on the lateral oscillatory stability of a 
hypothetical supersonic airplane. The investigation included an automatic 
pilot sensitive to a displacement in either yaw or. roll and an auto— 
matip pilot sensitive to either the yawing or rolling angular velocity. 

The calculations were made for each type of automatic pilot acting inde- 
pendently of the others. An idealized control system without lag was 
assumed for the calculations. 

The results of the investigation indicated that all the automatic 
pilots improved the stability of the original unstable Bhort— period oscil- 
lation. The only one of the automatic pilots which resulted in an oscil- 
lation that satisfied the NACA and military criterions for satisfactory 
damping— period relationship, however, is an automatic pilot sensitive to 
the yawing angular velocity and geared to the rudder so that rudder 
control is applied in proportion to the angular velocity. 


INTRODUCTION 


The lateral— stability boundaries calculated in reference 1 indicated 
that for high-speed airplanes designed with high wing loadings and srwept- 
back wings more directional stability is required for oscillatory 
stability than for an airplane with a straight wing and lower wing loading. 
Subsequent lateral stability investigations (references 2 and 3) showed 
that the directional stability required for oscillatory stability may 
be reduced if the principal longitudinal axis of the airplane is inclined 
above the flight path. Additional, unpublished, dynamic— stability calcu- 
lations on the effect of an automatic pilot which gives control pro- 
portional to the angular displacement in either yaw or roll or an automatic 
pilot which gives control proportional to either the yawing or rolling 
angular velocity alBo indicated that the use of the automatic pilot permits 
a reduction in the directional stability required for oscillatory stability. 
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These calculations for the effect of the automatic pilot were made at 
the time the theoretical investigation reported in reference 3 was carried 
out j hut the results of the calculations were not published then because 
the stabilizing effect that would be obtained by inclining the principal 
axis above the flight path was believed to be sufficient to obviate the 
necessity of installing an automatic pilot. Recent lateral— stability 
analyses of several airplanes have shown, however, that the oscillatory 
stability is not satisfactory even when the principal axiB is inclined 
above the flight path. The use of an automatic pilot therefore offers 
another means for improving the oscillatory stability. 
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SIMBOLS AND COEFFICIENTS 

angle of roll, radians 

angle of yaw, radians 

angle of sideslip, radians (v/V) 

yawing angular velocity, radians per second (d^/dt) 

rolling angular velocity, radians per second (d0/dt) 

sideslip velocity along the Y— axis, feet per second 

airspeed, feet per second 

mass density of air, slugs per cubic foot 

dynamic pressure, pounds per square foot 

wing span, feet 

wing area, square feet 

weight of airplane, pounds 

mass of airplane, slugs (w/ g) 

acceleration due to gravity, feet per second per second 
relative— density factor (m/pSb) 

inclination of principal longitudinal axis of airplane 
with respect to flight path, positive when principal 
axiB is above flight path at the nose, degrees 

angle of flight path to horizontal axis, positive in a 
climb, degrees 

radiuB of gyration in roll about principal longitudinal 
axis, feet 
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radius of gyration in yaw about principal -vertical axis, 
feet 


nondimensional radius of gyration in roll about principal 
longitudinal axis 

nondimensional radius of gyration in yaw about principal 
vertical axis ^b^/b^ 

nondimensional radius of gyration in roll about longi- 
tudinal stability axis / , — \ 

(n/%o 2cos2t 1 + K Zo 2sin2l l ) 

nondimensional radius of gyration in yaw about vertical 


stability axis 


0% o 2 cos 2 tj + Kx o 2 Bin 2 T] 


nondimensional product— of— inertia parameter 


((^Zo 2 ~ Klo 2 )sin n cos 
trim lift coefficient 


roUing-moment coefficient 


yawing-moment coefficient 


lateral— force coefficient 


Bolling moment 


'Lateral force 


effective— dihedral derivative , rate of change of rolling- 
moment coefficient with angle of sideBlin per radian 
(dCj/dP) 

directional— stability derivative, rate of change of 

yawing-moment coefficient with angle of sideslip, per 
radian (3Cn/8p) 

lateral— force derivative, rate of change of lateral— force 
coefficient with angle of sideslip, per radian 
(SCy/dP) 

damping— in— yaw derivative due to the airplane, rate of 
change of yawing-moment coefficient with yawing- 
angular— velocity factor, per radian ^^ n y^rb^ 
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damping— in— yaw derivative due to the automatic pilot 


rate of change of yawing-moment coefficient with rolling- 
angular— velocity factor, per radian (aory^) 

damping— in— roll derivative due to the airplane, rate of 
change of rolling-moment coefficient with rolling- 

angular— velocity factor, per radian 

damping— in-roll derivative due to the automatic pilot 

rate of change of lateral— force coefficient with rolling— 
angular-velocity factoij per radian 

rate of change of lateral— force coefficient with yawing— 
angular-velocity factoij per radian 

rate of change of rolling-moinant coefficient with yawing— 
angular-velocity factor, per radian 

rate of change of yawing-jnoment coefficient with angle 
of yavj per radian 



rate of change of rollin 

- „ fbQ 

of jav, per radian 



iment coefficient with angle 


rate of change of lateral— force coefficient with angle 

f *°A 

M ) 


of yawj per radian 



ment coefficient with angle 


rate of change of yawing 
' of roll, per radian 


rate of change of roUing-momant coefficient with angle 
of roll, per radian 
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rate of change of lateral-rforce coefficient with angle 
of roll, per radian ^ 

rate of change of yawing-moment coefficient with rudder 
deflection, per radian 


rate of change of rolling-^aom^nt coefficient with rudder 
deflection, per radian 


3/X 

te 

\^r> 


rate of change of lateral— force /coefficient with rudder 

© 


deflection, per radian I 


rate of change of yawing-moment coefficient with aileron 

[tea/ 


deflection, per radian 


rate of change of rollin 
deflection, per radian 



t coefficient with aileron 


rate of change of lateral— force coefficient with aileron 

/dOy\ 

deflection, per radian I J 


control-gearing ratio, rate of change of rudder deflec- 
tion with angle of yaw 

control-gearing ratio, rate of change of aileron deflec- 
tion with angle of roll 

control-gearing ratio, rate of change of rudder deflec- 
tion with yawing angular velocity 

control— gearing ratio, rate of change of aileron deflec- 
tion with rolling angular velocity 

thus, seconds 

nondimens ional time parameter based on span (Vt/b) 

differential operator 

d&b 
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El Eouth’s dlBcrlmlnant for a quintic equation 

E 2 Eouth’s discriminant for a quartic equation 

P period of oscillation , seconds 

Tl/2 time for amplitude of oscillation to change "by factor 

of 2 (positive value indicates a decrease to half 
amplitude., negative value indicates an increase to 
double amplitude) 

/g number of cycles required for amplitude of periodic 

mode to change by factor of 2 (positive value indicates 
a decrease to half amplitude; negative value indicates 
an increase to double amplitude) 

A,B.,C.,D.,E,,F coefficients of lateral— stability equations 


SCOPE OF INVESTIGATION 


The oscillatory— stability boundaries were calculated for a hypo- 
thetical airplane to show the effect on the oscillatory stability of an 
automatic pilot which gives control proportional to the angular displace- 
ment in either yaw or roll or an automatic pilot which gives control 
proportional to either the yawing or rolling angular velocity. The calcu- 
lations were made for each type of automatic pilot acting independently 
of the others. An idealized proportional control system in which lag 
effects were neglected was assumed for the calculations. The relation- 
ship of the calculated boundaries to the motion of the aircraft was 
investigated by determining the period and damping of the oscillatory 
mode and the damping of the aperiodic mode from the roots of the charac- 
teristic lateral— stability equation. 

The mass and aerodynamic parameters of the hypothetical airplane 
in the cruising condition are presented in table I. The results of the 
calculations presented are based on the assumption that the product of 
inertia is zero. 
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EQUATIONS OF MOTION 

The linearized equations of motion, referred to stability axes, 
for the condition of controls fixed are: 

Bolling 

+ ExzPb = c 2pP 

Yawing 

2^ K zV^ ^ZZpb 2 ^) = °n p P 

Sideslipping 

,2p.-b(l>bP + I>b t) = CypP 

+ ^ tan 7) y 

When 0 o e^ ^ is substituted for 0, i|f 0 e^ St for \}r, and p 0 e^ B ^ 
for P in the equations written in determinant form, \ muBt he a root 
of the stability equation 


+ + g c 2 r :D b 1 l f 

+ + i<V>bi 

+ TpXjPbtf + °L^ + ^Y^Tj; 


AX^ + BX 3 + cx 2 + DX + E = 0 

where 

A = 8^3 (k-j 2 ^ 2 - %z 2 ) 

B =-2^ 2 (2K x % 2 CY p + Kj 2 ^ + K Z 2c l p “ ^XZ^p “ K XZ G l r ~ %Z C n p ) 

C = Ub(Kx 2 Cn r CYp + ^^bEl 2 Qnp + K^CipCyp + |0n r C2 p - ExzC2 r CY p 
~ ^b%Z C 2p “ 2 C np C 2 r “ ^p^Z^p + %Z C n p C Y p “ E Z 2c Y p G 2p 

“ % 2c Y r °np + %Z°y r %) 
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D = - ^n r c Zp°Yp “ ^bCZpCnp + jQnp c 2 r c Tp + ^bQopCZp + 2 ^Ci£xzCn p 

- 2l^ b C L K: z S C Z p - 2^ b K x 2 C np C L tan 7 + tan 7 + 

- t%%°T r - JOvPnp^ + iqnpOipOjp 


E = iQL(Cla r O Jp - C 2 r Cn p ) + |cl tan r (%Cn f - % 0 lf~) 


If an automatic pilot sensitive to a displacement in yaw is 
installed in the airplane and rudder control is applied in proportion to 
the displacement, the stability derivatives Cj^, and Cy^ are 

introduced into equations (l) . These derivatives will he functions 
of Qng , C 2 g r j and Cyg^, respectively, and of the control-gearing 

8&r 

ratio - — : that is. 

8^ 


and 





_ 86 

C 7 — On — 

z 8 r 


85^ 


% ~ °*B r 


Normally for present-day airplanes, Cjg and C Y g are Bmall; 

equations (l) will therefore he simplified hy neglecting the terms 
involving and Gy^. If appears in equations (l), the 

stability equation becomes 

A\5 + BX 1 ^ + CX3 + DX 2 + EX + F = 0 


where the following terms are added to the terms of the coefficients 
presented previously: 


for 0 

- ^b 2 % °ny 
»*b C 2p9n + + 2 ^b K X 2c Y p G n t 


for D 
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for 1 


for F 


“ 2%%% + JCy^C^ 


If the automatic pilot is made sensitive to a displacement in roll 
and if aileron control is applied in proportion to the displacement, the 
derivatives Cj Cn^, and Cy^ are introduced into equations (1). These 

derivatives will he functions of C^g , Cng , and Cy 5 , respectively, 

' 85 

and of the control— gearing ratio — S. The derivative C^g and Cy~ 

op a 

were assumed to he small; therefore, the derivatives. Cn^ and Cy^ will 

he neglected and only the derivative will he taken into account in 

equations (1). The coefficients of the stability equation will then 
include the following additional terms: 


for C 


for D 


for E 




V 


+ 2 ^h K z 2 °yp c z^ 

~ | C n r C Yp C 2^ ~ 2 ^b °np C Z0 + 


for F 

Citan 7 CnpC^ 


For the cases of the automatic pilot designed to give control pro- 
portional to the yawing or rolling angular velocity, only the derivatives 
for damping in yaw and roll Cq^ and were assumed to he effectively 

increased. Actually, the derivatives C , Cy r , C^, and Cy^ are 

affected hy this type of automatic pilot. . The effect of the automatic 
pilot on C Cy r , Qnp, and Cy^, however, was neglected because these 

derivatives are functions of C|g , Cy^ , ^S a J reB ^ ec ^ ve ^-yj 

nil of which are nor mall y smal 1 for present-day airplanes!, The deriva- 
tive C was effectively increased when the rudder was assumed deflected 
^r 
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in proportion to the yawing angular velocity , whereas C 7 was effect— 

4 p 


ively increased when the aileron 
to the rolling angular velocity. 

due to the automatic pilot are 


was assumed deflected in proportion 
The expressions for and 


and 


ACn r - Qn^ 


8&r 2V 



^a 2Y 

80 "b 


EBSDLTS AND DISCUSSION 


.Figure 1 clearly shows that the probable range of Cy^, G\^ combi- 
nations for the hypothetical high-speed airplane is located almost 
entirely in the oscillatorily unstable region. This boundary represents 
the case in which the principal longitudinal axis iB alined with the 
flight path; therefore, the product of inertia is zero. The effect of 
automatic stabilization on the neutral-oscillatory— stability boundary iB 
presented in figures 2 to 5 . In each figure the ordinate is the 
directional— stability derivative Cnp and the abscissa is the stability 
derivative introduced by the automatic pilot. The curveB on each figure 
represent three different values of the effective-dihedral derivative Cjp. 
The control-gearing ratio is also plotted as abscissa in figures 2 to 5 
to indicate the relation between the control-gearing ratio and the 
stability derivative introduced by the automatic pilot. The control- 


gearing ratios 


that 


and 


8S r 85 a ctfl,, 
8i{r 80 8ijr 


and 


85, 

80 


were calculated on the assumption 


°BS r “ ~ 0 ' 1 


' j 5 £ 


= - 0.1 


Automatic Pilot Sensitive to Displacement in Yaw 

If the automatic pilot is sensitive to a displacement in yaw and 
rudder control is applied in proportion to the displacement, the deriva- 
tive Cn^ is introduced into the equations of motion. The lateral- 

stability equation now becomes a quintic equation as shown in the previous 
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section entitled "Equations of Motion". The necessary and sufficient 
conditions for complete stability for a system which resultB in a quintic 
stability equation are derived in reference 4 and are also presented in 
reference 5* A simple derivation of the necessary and sufficient condi- 
tions for neutral oscillatory stability and for complete stability is 
given in the appendix. The necessary and sufficient conditions for 
neutral oscillatory stability are that 

Ep = (BC - AD)(DE - CP) - (BE - AF) 2 = 0 


and that BE — AF and BC — AD are of the same sign. The necessary 
and sufficient conditions for complete stability are that the coeffi- 
cients A, B, D, and F be positive, Bp > 0* a* 1 * 3 - BE — AF > 0. 


Figures 2(a) to 2(c) show the neutral-oscillatory-etability boundary 


Bp = 0 and the curve BE — AF = 0 plotted as a function of and 

On the shaded 


*1 

or 






for C lp = -0.10, -0.30, and -0.50, respectlTel,. 


side of BE — AF = 0, the airplane has at least one mode of motion which 
is unstable. The significance of the boundary El = 0 located on the 
shaded side of BE — AF = 0 iB still valid; that is, a stable oscilla- 
tion becomes unstable upon crossing the boundary. A clearer relationship 
between the boundary Ei = 0 and the stability of the airplane motion 
may be obtained from a study of the damping and the period of the oscilla- 
tion and the damping of the aperiodic modes obtained from the roots of 
the quintic lateral-stability equation. The reBultB are presented in 
table H(a) for the case of = -0.1 at C n p = 0.15 and C n p = 0-55 


for several values of Cq^. p or point A in figure 2(a), which corre- 
sponds to the point at the values of C = -0.1, C n p = 0.15, 

and = 0, the roots of the stability equation indicate an unstable 

short-period oscillation, two subsiding aperiodic modes, and a zero root. 
The zero root means that the airplane iB insensitive to displacements in 
yaw. As is increased negatively from point A to point B, a stable 

long— period oscillation is introduced in addition to the unstable short- 
period oscillation and one subsiding aperiodic mode. Passing through 
the boundary Ep = 0, from point B to point C, causes the stable long- 
period oscillation to become unstable, but does not appreciably affect 
the other modeB of motion. Two unstable oscillatory modes of motion now 
exist. As is further increased negatively to point D, the insta- 

bility of the short— period oscillation is reduced; whereas the instability 
of the lo ng— period oscillation increases and its period is shortened. 
Crossing th rough the boundary Bp - 0 to point E causes the short— period 
oscillation to become stable; however, as continues to increase 

negatively, the damping of the short-period oscillation does not improve 
sufficiently to meet criterion^ for satisfactory dynamic stability (as 
described subsequently). This conclusion for large negative values 
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of where cross— coupling effects are negligible, could be checked by 

considering C n ^ — Cjj to be the spring constant in the mass-spring dashpot 
system as expressed by the one— degrea-of— freedom equation in yaw 


2 Mb^Z 2;D b 2 ~ |Cn r »b ~ Cn p )]t = 0 


An increase in the constant term Cn^ - Cn p of the equation reduces the 
period of the oscillation, but does not affect the damping. 


The results presented in table H(a) for C 2 p = -0.1 and Qnp = 0.55 

when used in conjunction with figure 2 (a) indicate that upon crossing 
through the boundary Hi = 0 the long— period oscillation becomes unstable. 
As Qmjr i 0 increased negatively, both the damping and the period of the 
short— period oscillation vary slightly and in such manner that the number 
of cycles required to damp to half amplitude Cw 2 is approximately 


constant. With increasingly negative values for Cn^-, the long— period 
oscillation becomes more unstable and its period is reduced. 


The results of these calculations for the hypothetical airplane, 
therefore, indicate that an automatic pilot sensitive to a displacement 
in yaw does not substantially improve the stability of the Bhort -period 
oscillation. 


Automatic Pilot Sensitive to Displacement . in Eoll 


If the automatic pilot is sensitive to a displacement in roll and 
aileron control is applied in proportion to the displacement, the 
derivative is introduced into the equations of motion. In the 

present analysis, level flight is assumed or 7 = 0° and the stability 
equation is, therefore, a quartic. The necessary and sufficient condi- 
tions for neutral oscillatory stability are that 


E 2 = BCD - AD 2 _ B 2 E = 0 


and that the coefficients B and D are of the same sign (reference 6 ). 


The neutral— oscillatory— stability boundary Bp = 0 is presented 
in figure 3 for = -0.10, -0.30, and -0.50. The complete curves 

of E 2 = 0 for C ^ = -0.30 and -O. 5 O, however, are not shown in 
figure 3 . Actually, the complete curves for these values of C 7 ,, are 




similar in shape to the curve for 


G Z p = -0-10 


and intersect the line 
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of 


= 0 at values of greater than those plotted in figure 3. 


The houndaries B2 = 0 in figure 3 indicate that when the oscillation 
is unstable for = 0 , oscillatory stability is obtained as 

exceeds the value corresponding to a point on the boundary Eg = 0 * 

If the oscillation is stable at = 0 , oscillatory instability will 

occur for a limited range of values of C as prescribed by the 
boundary E2 = 0 . Further negative increases in C^, however , will again 


result in oscillatory stability. The damping and period relationship of 
the oscillatory modes for C^p =-0.10 and Qn.g = 0.15 and 0.^5 are 

presented in table H(b). For point A in figure 3, which corresponds to 
the point at the values of C n p = 0.15 and C = 0, the roots of the 

stability equation indicate an unstable short— period oscillation and two 
subsiding aperiodic modes. As Cirf is increased negatively to point B, 


the two subsiding modes combine to form a stable short— period oscillation. 
Between points A and B, the boundary for two equal roots occutb, beyond 
which two oscillations exiBt (reference 6). Upon passing through thiB 
boundary of equal roots, the period of the newly formed oscillation is 
very long, but it rapidly decreases as indicated by the resultB of the 
calculations at point B where the period iB approximately 4 BecondB. 
Crossing through the boundary B2 = 0 , from point B to point C, causes 


the original oscillation to become stable. Ab C z ^ continues to increase 

negatively, the period of one of the oscillations remains constant and 
its damping decreases slightly whereaB the period of the other oscillation 
decreases and its damping increases slightly. For both oscillations, the 
number of cycles required to damp to half amplitude C increases 

as increases. For large negative values of C^, the period and 

damping of each oscillation may be approximated by the one-degree-of— 
freedom equations of motion in roll and yawj 


Bolling 

^(ibEz 2 !!) 2 - |CZpBb - - 0 

Yawing 

(fyftfab 2 _ Icn^ _ = 0 

where 





Thus it is seen that as C ^ increases negatively, the system acts as 

if two independent mass— spring daBhpot systems were in operation without ■ 
any cross-coupling effects. The derivative acts in the capacity 
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of a spring constant and tlie period, therefore, varies inversely with Cjj. 
A comparison of the results in table H(b) for Cn^ = O.I 5 and Cnp = 0.45 
at C = — 0 . 0 5 , —0.10, and -0.20 clearly shows that the oscillation 
described by the rolling equation is approximately independent of Cnp • 

For small negative values of ^ 1^3 the equations of motion for the three 

degrees of freedom involved in lateral motion must be solved simultane- 
ously to determine the period and damping of the modes of motion. 

Thus an automatic pilot sensitive to an angular displacement in roll 
stabilizes the original unstable oscillation and introduces an additional 
stable oscillatory mode with a shorter period than the original oscil— ' 
lation; however, throughout the range of C20, the period and damping 

relationship of at least one of these oscillations may be objectionable 
to the pilot. 


Automatic Pilot Sensitive to Rate of Displacement 

Another type of automatic pilot included in the investigation was 
one sensitive to either the yawing or the rolling angular velocity. Rudder 
control was assumed to be applied in proportion to the yawing angular 
velocity, thereby increasing the damping in yaw derivative C nr ; whereas, 

aileron control was assumed to be applied in proportion to the rolling 
angular velocity causing an increase in the damping— in— roll derivative c 2p. 
For both cases the stability equation is a quartic and the necessary and 
sufficient conditions for neutral oscillatory stability are similar to 
those conditions described in the previous discussion of an automatic pilot 
sensitive to roll. 

Automatic pilot sensitive to yawing angular velocity.— Figure 4 shows 
the neutral— oscillatory-stability boundaries R 2 = 0, for C}p = -0.10, 

— 0 . 30 , and — 0 . 50 . The boundaries indicate that for certain values of 

oscillatory stability is obtained provided a definite lower limit of ACn r 
iB exceeded. However, as Z£n r is further increased, a critical value 

is reached beyond which the airplane is unstable. This result is caused 
by the fact that if the damping in yaw is made sufficiently large, the 
airplane motion in yaw is restricted. If the motion is then analyzed on 
the assumption that only two degrees of freedom (roll and sideslip) remain, 
oscillatory instability will occur for negative values of (positive 

effective dihedral). (See reference 3 .) P 

The period and damping relationship for = — 0.10 and C n p = 0.15 

for several values of £&a r are presented in table n(c) . At point A in 

figure 4, which corresponds to the airplane without an automatic pilot, a 
short— period unstable oscillation and two subsiding aperiodic modes are 
obtained. When the derivative due to the automatic pilot AC nr is intro- 
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duced or a shift from point A to point B occurs, the instability of the 
short-period oscillation is improved. Crossing through the boun dar y Bg = 0 

to point C causes the unstable short— period oscillation to become stable. 

A heavily damped long-period oscillation which is formed from the combi- 
nation of the two subsiding aperiodic modes also appears for conditions 
corresponding to point C. As ACn r is increased negatively between 
point C and point D, the period of the short— period oscillation increases 
and the damping is improved.; thereby greatly reducing the number of cycles 
required to damp to half amplitude as expressed by the value of C-jyg. 

The effect of increasing AC ^ negatively on the long— period oscillation 

is to reduce the period and to decrease the damping. This oscillation 
with a period of about 7*5 Beconds becomes unstable upon passing through 
the boundary Eg = 0 to point D. Tor ACn r = - 44 . 0 , the rootB of the 

stability equation show that the stable oscillation breakB down into two 
heavily damped subsiding aperiodic modes and the unstable oscillation 
becomes more unstable. Because of the comparatively long period of the 
unstable oscillation, however, pilots might not find thiB type of insta- 
bility difficult to control. 


Thus the results indicate that an automatic pilot sensitive to the 
yawing angular velocity and geared to the rudder causes a marked improve- 
ment in the original unstable short-period oscillation without introducing 
any additional modes of motion which might be objectionable to the pilot. 


Automatic pilot sensitive to rolling angular velocity .— The neutral— 
oscillatory-stability boundary I?2 = 0 for the airplane equipped with 


an automatic pilot sensitive to rolling angular velocity is presented in 
figure 5 f°n several values of C n . The boundaries indicate that oscil— 

latory stability is obtained provided a definite lower limit 'of AC7 Is 

fc p 

exceeded. An ex am i na tion of the period and damping of the oscillation, 
presented in table H(d) for = —0.1 clearly shows that, although 

the unstable short— period oscillation does become stable as ACjp 1 b 

increased negatively, any additional damping in roll introduced into the 
system does not improve the damping of the oscillation sufficiently to 
result in a satisfactory damping and period relationship. The effect of 
adding more AC 2 into the system is simply to increase the damping of 


one of the subsiding aperiodic modes while causing a reduction in the 
damping of the other subsiding aperiodic mode. It is important to note 
that for some airplane configurations adding more damping in roll into 
the system will cause the oscillatory mode to become less stable. 
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Effect of Automatic Stabilization on the CriterionB for 
Satisfactory Damping-Period Belationship. 
of the Oscillatory Mode 

The curves that define the satisfactory damping-period relationship 
of the oscillatory mode as a function of the period are presented in 
figure 6. The dashed curve represents the NACA criterion (reference 7 ) 
and the solid curve represents the military criterion (references 8 and 9) . 
Both curves agree very well with each other for periods greater than 
6 seconds whereas for periods from 0.8 to 6 Beconds, the military criterion 
is more stringent than the NACA criterion. For periods smaller t hen 0.8 
of a second, the NACA criterion is more conservative than the military 
criterion. 


The relative merits of the several types of automatic pilots dis- 
cussed in the present paper can he clearly seen by comparing the damping- 
period relationship of the oscillatory mode, as affected by the particular 
type of automatic pilot, with the criterions shown in figure 6. Curves 

of — - — are presented as a function of the derivative introduced by the 
C l/2 


automatic pilot for an automatic pilot sensitive to either a displacement 
in yaw or roll in figure 7 and for an automatic pilot sensitive to either 

the yawing or rolling angular velocity in figure 8. One curve of ~ - 

c l/2 


is presented for each of the automatic pilotB investigated despite the 
fact that two oscillatory modes may occur, Bince the airplane motion would 
'be objectionable to the pilot if only one of the oscillatory modes does not 
satisfy the criterion for the damping-period relationship. An improvement 
in the original damping— period relationship as each particular type of 
automatic pilot is introduced is noted in figures 7 and 8 by the fact 

that n =7 changes fhom. a negative to a positive sign. For an automatic 
^ 1/2 

pilot sensitive to a displacement in yaw (the solid curve in fig. 7) the 
values of correspond to point b located in the unsatisfactory 

region of figure 6 according to both criterions. For an automatic 
pilot sensitive to a displacement in roll (the dashed curve in fig. 7 ) 
the oscillatory mode satisfies the NACA criterion for values of C 


from - 0.015 to - 0 . 05 . However, for the period of the oscillation during 
this interval, which ranges from a period of 1.4 to about 2.5 seconds, the 
military criterion is not satisfied. The values of ^ for the dashed 


curve in figure 8, which corresponds to an automatic pilot sensitive to 
the rolling ang ular velocity, indicate that neither one of the criterions 
in figure 6 is satisfied. The solid curve in figure 8, which corresponds 
to an automatic pilot BenBitive to the yawing angular velocity, shows 
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that as AC^ Increases negatively, the oscillatory mode satisfies the 

NACA and military criterions. An additional oscillatory mode exists 
which also satisfies the criterions for the range of AC^ shown in the 

figure. However, as ACn r continues to increase negatively to a value 
of approximately -25*0, the value of C-^ becomes negative. (See 
table 3Z(c) . ) 

Figures 7 and 8 indicate that for thiB hypothetical supersonic 
airplane an automatic pilot sensitive to the yawing angular velocity, 
which effectively increases the derivative C n ^, is the most desirable 

type of automatic pilot to be used to obtain a satisfactory damping — 
period relationship of the oscillatory mode. 


CONCLUSIONS 


The following conclusions were drawn from a theoretical investi- 
gation carried out to determine the effect of automatic stabilization on 
the lateral oscillatory stability of a hypothetical supersonic aircraft: 

1. An automatic pilot sensitive to a displacement in either yaw or 
roll and an’ automatic pilot sensitive to either the yawing or rolling 
angular velocity improved the damping-period relationship of the original 
unstable short-period oscillation. 

2. The only one of the several types of automatic pilots investi- 
gated which resulted in an oscillation that satisfied the NACA and 
military criterions for satisfactory damping— period relationship is an 
automatic pilot sensitive to the yawing angular velocity and geared to 
the rudder so that rudder control is applied in proportion to the angular 
velocity. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya. , December 9, 19^8 
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APPENDIX 

DERIVATION OF THE NECESSARY AND SUFFICIENT CONDITIONS FOR 
NEUTRAL OSCILLATORY STABILITY AND COMPLETE STABILITY 
■ OF THE QUINTIC EQUATION 
By Leonard Sterafield and Ordway B. Gates, Jr. 


The necessary and sufficient conditions for neutral oscillatory 
stahility of the quintic equation are that the coefficients of the 
stability equation 

A\5 + BX^ + CX 3 + DX 2 + EX + F = 0 (Al) 

satisfy Routhls discr iminant set equal to zero 

R-L = (BC - AD )(DE - CF) - (BE - AF) 2 = 0 

and that BC — AD and BE — AF have the same Bign. The expression 
for R^ = 0 can be derived by assuming that the quintic equation has two 
roots X = ±icn, where to is the angular frequency of the neutrally stable 
oscillation. This assumption is based on the fact that for the condition 
of neutral oscillatory stability the real part of the complex root muBt 
be zero. If X = loo Is substituted in the equation (Al) the following 
two expressions are obtained: 

Ao5 — 0 cd3 + Ed = 0 (A2) 

Bcd 1- - Dm 2 + F = 0 (A3) 


Solving equations (A2) and (A3) simultaneously, thereby eliminating the 
aft terms, gives the expression 


CD 2 _ BE - AF 
BC - AD 


(A4) 


Substituting equation (Ak) into either equation (A2) or (A3) results in 
Routh* b discriminant 

(BC— AD) (DE— CF) - (BE-AF) 2 = 0 


It is seen from equation {Ah) that the expression 


_ BE - AF 
NBC - AD 


co 
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defines the angular frequency of the neutrally stable oscillation. The 
symbol cd represents the frequency of the neutrally stable oscillation 
only if BC — AD and BE — AE are of the same sign, since cd must have 
a real value if the root iB to represent an oscillation. If BC — AD 
and BE — A3? are of different sign and Ei = 0 is satisfied, cd 1b an 
Imaginary quantity and the two roots of the quintic equation given 
by X = ±icD are two real roots equal in magnitude but opposite in Bign. 

The necessary and sufficient conditions for complete stability of 
the quintic equation are derived in reference h and are presented in a 
condensed form in reference 5* A much simpler method for obtaining 
these same conditions is presented in the following analysis. 

Assume that the quintic equation 

AxP + BX^ + CX^ + DX*~ + EX + F — 0 


has the roots 


Xq — aj_ + ibq 

X2 = aj — ibq 

X 3 = ag + ibg 
X||_ = — 11*2 

X 5 = a 3 


(A5) 


Between the coefficients and roots of the quintic equation the following 
relationships exist: 


1 

(A6a) 

— 2 (aq + a g ) — a 3 

(A6b) 

aq2 + 13^2 + ^2 + tg 2 + 4a]_a2 + 2a 3 ^aq + 

(A6c) 

— S^a^^ag^ + b2^^ + ^aq^ + bq^J 


— a 3 (aq 2 + aQ 2 + bq 2 + b2 2 + 

(A6d) 

(ai 2 + bi 2 )(a 2 2 + b 2 2 ) 


+ 2a 3 jaq (a^ 2 + b£ 2 ) + a^aq 2 + bq 2 )] 

(A6e) 
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F = -a 3 (a 2 2 + 'bi 2 )(a 2 2 



(A£f) 


The requirement for complete stability is that a-^, ag, and a^ 

be negative. Substitution of negative values of a^, a^ 3 and makes 

all the coefficients positive; therefore, the first condition for 

complete stability is that all coefficients be positive. From equation (AfSf) 

it can be seen that the F coefficient will be positive only when a^ is 

negative. This first condition, therefore, only insures the stability 
of the root = a^. Additional requirements are needed to make a^ 

and ag negative. On the boundary between stability and instability, 

where = 0, the equations ( A6 ) become 

A 
B 

C 
D 
E 
F 

that is, five equations in the four variables a^, a^, b^, and bg. 
Elimination of these variables leads to Routh ? s discriminant 

Rq = (BC - AD)(DE - CF) - (BE - AF) 2 = 0 


1 

- 2a l“ a 3 

a^ 2 + b^ 2 + bg 2 + 2a„a. 


'3 1 


-2a 1 b 2 2 - a 3 (a 1 2 + b^ + bg 2 ) 
( a l 2 + t l 2 ) t 2 2 + 2a 3 aib 2 2 


= — a 0 ( a ! 2 + ^ l 2 )^ 2 


*3 VI 


(A7) 


which was derived by a much simpler procedure in the preceding analysis 
of the conditions necessary for neutral oscillatory stability. It is 
apparent that this discriminant can be equal to zero only if BC — AD 
and DE — CF are of the same sign. When ag = 0 the following relation- 
ship can be obtained from equations (A7): 

BC - AD = -2aq^a 3 2 + aq 2 + bq 2 ) - ka^aq 2 


Since a^ must be negative if all coefficients are positive, BC - AD 
must be positive if a^ is negative. Therefore, DE — CF must also be 
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positive if the discriminant is to "be equal to zero, which implies 

that BC > AD and DE > CF. Then C > ^ and c < BE, also M>AD. 

B- F F B 1 


therefore , 


BE - AF > 0 


(AS) 


Thus, for complete stability a second necessary condition is that BE — AF 
must he greater than zero. 


It can he shown from the simultaneous solution of equations (A7), 
where = 0, that 

BE - AF = -fia-L 



Since the terms within the bracket are all positive, BE — AF > 0 if, 
and only if, a-^ is negative. This verifies the correctness of 

equation (a8) . 

In order to determine the sign of Routh’B discriminant for a condi- 
tion of complete stability, it is only necessary to substitute the roots 
of a completely stable case into equations (A 6) and form the discriminant. 
This substitution establishes the criterion that for complete stability 

% = (BC - AD)(DE - CF) - (BE - AF) 2 >0 (A9) 

In summary, the necessary and sufficient conditions for complete 
stability of the quintic equation are that all coefficients muBt be 
positive, BE — AF > 0, and Routh's discriminant muBt be positive. The 
conditions obtained as a result of this analysis agree with the conditions 
as presented in reference 5 . The conditions Btated in reference 5 ar© 
entirely adequate although they do not specifically state that the coeffi- 
cients C and E must be positive, since the condition that BE — AF 
must be greater than zero inherently demandB that these two coefficients 
"be positive as is evident from the derivation of this condition in the 
preceding analysis. 
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TABLE I.- STABILITY DERIVATIVES AND MASS CHARACTERISTICS 
OF HYPOTHETICAL AIRPLANE IN CRUISING CONDITION 


v/s, lb /ft 2 ....... 

"b ,ft 

p 3 slugs/ft^ . 

Y 3 ft/aec ........ 

7, deg 

QL ........... 

........... 

kjo, ft ........ . 

ft 

deg 

0i VJ per radian . . . . . 
C^^, per radian 

Cq , per radian ..... 
P 

C^, per radian 

Cyp, per radian 

Oy^, per radian . . . . . 
Cyp, per radian 

S(fuaelage)' ® er radlan 

Cnpj per radian . . . . . 

S(tail)' Ver radian . . 
C^pj per radian . . . . . 


»..«.*«. ^1-0 

20 

0.0002 

O O O . . O . 1^-65 

0 

...... 0.372 

620 

....... 2.02 

9‘Sb 

0 

...... -0.197 

0.0929 

..... -0.00732 

• ~ lA7 ^(tail) 

. 0 

0 

• _1 ' 33 S(tall) 
...... -0.25 

Variable 

^P ^Hp ( fus elage ) 
Variable 
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TABLE II.- PERIOD AND DAMPING OF THE OSCILLATORY MODES 
AND DAMPING OF THE APERIODIC MODES 


[c, = -O.l” 

L 2 P J 


(a) Automatic Pilot Sensitive to Displacement in Yaw- 


Points Cnp Cn ^ 



Oscillatory mode 


p T l/2 Cl/2 


- 2.11 


0.15 0 


.15 -. 



.15 -.0035 
.15 -.13 
.15 -.15 


.15 -. 


r 3.61 -7.88 
L59.2 247.0 
r 3.61 -7.97 
144.2 -835.O 
j 2.73 -217.0 


110.3 
‘ 2.62 
10.0 
r 2.45 
■ 9.^0 
J1.96 
18.65 
1.95 


-3.75 

330.0 

-3.53 

68.0 

-3.22 

26.7 

-2.72 

11.6 


.55 0 


-.0035 


1.95 

188.5 

' 1.95 
' 66.8 
1.93 

'29.7 

r 1.79 
' 15.6 
‘ 1.70 

12.5 

' 1.50 
.10.5 


11.6 

346.0 

11.6 

-825.0 

11.2 

-26.5 

10.2 

-7.28 

9-7 

-4.9^ 

9.0 

-3.71 



5-95 

3.92 

5.95 

-12.4 

5.80 

-.892 

5.70 

-.466 

5.70 

-.395 

6.0 

-•353 


0.827 

32.7 


1.06 

58.3 
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TABLE TT.— Continued 


(c) Automatic Pilot Sensitive to Yawing Angular Telocity 


Point b 

% 

ACn r 

88 r 

8* 

Oscillatory mode 

Aperiodic 

mode 

T l/2 

P 

t 1/2 

C l/2 






3.62 

-7.65 

-2.11 

0.827 

A 

0.15 

0 

0 

« 




32.70 






3.66 

-14.96 

-4.09 

.863 

B 

.15 

-•733 

• 05 





11.60 





f 

3-84 

4.38 

1.14 


C 

.15 

-4.40 

.30 

{: 

L79-90 

1.51 

.0084 







[4.13 

2.11 

.511 



• 15 

-7.33 

.50 


I16.OO 

1.45 

.0906 







J5.68 

.647 

.114 



.15 

-14.70 

1.0 


17.43 

3.44 

.463 







[6.69 

.432 

.0646 



.15 

-20.50 

1.4 


’ 7-42 

17.20 

2.32 







19-54 

.356 

.0373 


D 

• 15 

-25.40 

1.73 


'.7.^5 

-33-20 

-4.47 






r 




.483 


.15 

-44.0 

3-0 

j 

7.49 

-4.95 

-.661 

.135 
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TABLE II.— Concluded 


(d) Automatic Pilot Sensitive to Polling Angular Telocity 


Cnp 

A % 

86 a 

Oscillatory mode 

Aperiodic 

mode 

t 1/2 

P 

T l/2 

c l/2 




I 

'3-62 

-7.65 

-2.11 

0.1 7 

0.15 

0 

0 

1 




32.70 





‘3-59 

-20.90 

-5.82 

.345 

.15 

-0.44 

•03 





79.70 





r 3-63 

39-40 

10.90 

.172 

• 15 

-1.17 

.07 





157.50 




J 

tt.65 

26.50 

7.26 

.142 


-1.47 

.10 

< 

L 



189.40 





■3.71 

14.70 

3.96 

.052 

• 15 

-4.40 

.30 





501.60 





'3.72 

13.90 

3.74 

.032 

• 15 

-7.33 

.50 


— 



815.50 





3.74 

13-40 

3.58 

.016 

.15 

-14.70 

1.0 


. 



1593.0 















on the lateral oscillatory stability 
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